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The conventional order parameters in quantum matters are often characterized by 
'spontaneous' broken symmetries. However, sometimes the broken symmetries may 
blend with the invariant symmetries to lead to mysterious emergent phases. The heavy 
fermion metal URu2Si2 is one such example, where the order parameter responsible 
for a second-order phase transition at 7^=17.5 K has remained a long-standing mys- 
tery. Here we propose via ab-initio calculation and effective model that a novel spin- 
orbit density wave in the /-states is responsible for the hidden-order phase in URu2Si2. 
The staggered spin-orbit order 'spontaneous' breaks rotational, and translational sym- 
metries while time-reversal symmetry remains intact. Thus it is immune to pressure, 
but can be destroyed by magnetic field even at T = K, that means at a quantum 
critical point. We compute topological index of the order parameter to show that the 
hidden order is topologically invariant. Finally, some verifiable predictions are pre- 
sented. 
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Most states or phases of matter can be described by local order parameters and the associated 
broken symmetries in the spin, charge, orbital or momentum channel. However, recent discoveries 
of quantum Hall states, [3] and topological insulators (21 have revamped this conventional view. 
It has been realized[2 -4] that systems with combined time-reversal (T1Z) symmetry and large spin- 
orbit (SO) coupling can host new states of matter which are distinguished by topological quantum 
numbers of the bulk band structure rather than spontaneously broken symmetries. Subsequently, 
more such distinct phases have been proposed in the family of topological Mott insulators, [5 J 
topological Kondo insulators, [6 j topological antiferromagnetic insulators, [7 j. In the latter cases, 
the combined many-body physics and T1Z symmetry governs topologically protected quantum 
phases. Encouraged by these breakthrough developments, we search for analogous exotic phases 
in the heavy fermion metal URu 2 Si 2 , whose low-energy / states accommodate T1Z and strong 
SO coupling. This compound also naturally hosts diverse quantum mechanical phases including 
Kondo physics, large moment antiferromagnetism (LMAF), mysterious 'hidden-order' (HO) state, 
and superconductivity. H]] . 

In URu 2 Si 2 the screening of /-electrons due to the Kondo effect begins at relatively high 
temperatures, ushering the system into a heavy fermion metal at low-temperature. B Below 
7^=17.5 K, it enters into the HO state via a second-order phase transition characterized by sharp 
discontinuities in numerous bulk properties. lfTTT - [T3l [T6l The accompanying gap is opened both 
in the electronic structure ll9l [T4l - [T6l as well as in the magnetic excitation spectrum,[17j suggest- 
ing the formation of an itinerant magnetic order at this temperature. However, the associated 
tiny moment (~ 0.03/x#) cannot account for the large (about 24%) entropy release ifTOll and other 
sharp thermodynamic HH [16]] and transport anomalies [Q21 [13]] during the transition. Further- 
more, very different evolutions of the HO parameter and the magnetic moment as a function of 
both magnetic fieldflTH [201 and pressure[]2ll [22]] rule out a possible magnetic origin of the HO 
phase in this system. Any compelling evidence for other charge, orbital or structural ordering 
has also not been obtained. [23] Existing theories include multiple spin correlator [24 j, Jahn-Teller 
distortions [[25]], unconventional spin-density wave,(26l[27]], antiferromagnetic fluctuation, (28]] or- 
bital order, EUll . helicity order, [[29]], staggered quadrupole moment [7 j, octupolar moment, [31] hex- 
adecapolar order, [5 j, linear antiferromagnetic order. Il33ll incommensurate hybridization wave, [6] 
spin nematic order, [35] modulated spin liquid. Il36ll j-j fluctuations, (37 ], unscreened Anderson 
lattice model. Il38ll among others [8]. However, a general consensus for the microscopic origin of 
the HO parameter has not yet been attained. 



3 



Formulating the correct model for the HO state requires the knowledge of the broken sym- 
metries and the associated electronic degrees of freedom that are active during this transition. A 
recent torque measurement on high quality single crystal sample reveals that the four-fold rota- 
tional symmetry of the crystal becomes spontaneously broken ll23ll at the onset of the HO state. 
Furthermore, several momentum-resolved spectroscopic data unambiguously indicate the pres- 
ence of a translational symmetry breaking at a longitudinal incommensurate wavevector Q h = 
(1 ±0.4, 0, 0).[|T0l[T2l[T4l[T6J [Previous first-principle calculation has demonstrated that an accom- 
panying commensurate wavevector Q 2 = (1, 0, 0) might be responsible for the LMAF phase. Il33ll 
which is separated from the HO state via a first order phase transition. Il8l [T9l - [22l . As it is often 
unlike to have two phases of same broken symmetry but separated by a phase boundary, we expect 
that LMAF and HO phases are different.] In general, the order parameter that emerges due to a 
broken symmetry relies incipiently on the good quantum number and symmetry properties of the 
'parent' or non-interacting Hamiltonian. In case of URu 2 Si 2 , spin and orbital are not the good 
quantum numbers, rather the presence of the SO coupling renders the total angular momentum to 
become the good quantum number. Therefore, SU (2) symmetry can not be defined for spin or or- 
bital alone, and the 'parent' Hamiltonian has to be defined in SU (2) ® SU (2) representation. The 
'parent' Hamiltonian also accommodate other symmetries coming from its crystal, wavefunction 
properties which we desire to incorporate to formulate the HO parameter. 



RESULTS 



Ab-initio band structure. In order to find out the symmetry properties of the low-lying states, 
we begin with investigating the ab-initio 'parent' band dispersion and the FS of URu 2 Si 2 [|40l|42 
in Fig. [TJ The electronic structure in the vicinity of E F (±0.2eV) is dominated by the 5/ states 
of U atom in the entire Brillouin zone. (El [T414T61 l33l l42ll Owing to the SO coupling and the 
tetragonal symmetry, the 5/ states split into the octet J-\ (T 8 ) states and the sextet J=| (T 6 ) 
states. 03]] URu 2 Si 2 follows a typical band progression in which the T 8 bands are pushed upward 
to the empty states while the T 6 states drop to the vicinity of E F . The corresponding FS in Fig.[T]i 
reveals that an even number of anti-crossing features occurs precisely at the intersection between 
two oppositely dispersing conducting sheets. Unlike in topological insulators, 014]] the departure 
of the band crossing points from the TT^-invariant momenta here precludes the opening of an 
inverted band gap at the crossings, [2] and Dirac-cones crop up with Kramer's degeneracy in the 
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bulk states. Therefore, URu 2 Si 2 is an intrinsically trivial topological metal above the HO transition 
temperature. 

The SO interaction introduces two prominent FS instabilities at Q 2 = (1, 0, 0) and at Q h = 
(1 ± 0.4, 0, 0). The commensurate wavevector Q 2 occurs between same orbitals. Therefore, if 
this instability induces a gap opening, it has to be in the spin-channel, which is prohibited by T7Z 
symmetry and strong SO coupling. We argue (see Supplementary Information (SI) for details), in 
accordance with an earlier calculation, [33] that this instability is responsible for the LMAF phase. 
On the other hand, the incommensurate one, Q h , occurs between two different orbitals, and can 
open a gap if a symmetry between these orbitals and spins are spontaneously broken together. In 
other word, since spin-orbit coupling is strong in this system, individual spin- or orbital-orderings 
are unlikely to form unless interaction can overcome the spin-orbit coupling strength. On the other 
hand, a SO entangled order parameter in the two-particle channel can collectively propagate with 
alternating sign in the total angular momentum at the wavelength determined by the modulation 
vector. This is the guiding instability that drives spontaneous rotational symmetry breaking, while 
the T1Z symmetry remains intact (see Fig. [2^). This is because, both SU (2) groups for spin and 
orbital separately are odd under TTZ 9 but their product SU (2) ® SU (2) becomes even. As the 
parent state is not a non-trivial topological phase, a gap is opened to lift the FS instability. 

Low-energy effective model. Motivated by the above-mentioned experimental results and 
band structure symmetry properties, we formulate a simple and unified model by using the theory 
of invariants [44 J. We restrict our discussion to the low-lying T 6 bands and neglect the unfilled T 8 
bands. Due to j-j SO coupling and T1Z symmetry, the T 6 atomic states consist of three doublets, 
characterized by up and down 'pseudospins': mj=±|, ±|, ±|, where mj is the z component 
of J. On entering into the HO state, the FS instability commences in between the two doubly 
degenerate |mj|=§ and \ states only. 031 l35ll If no other symmetry is broken, the degenerate 
\mj\-\ state remains unaltered in the HO state, [44 j and hence they are not considered in our 
model Hamiltonian. Throughout this paper, we consistently use two indices: orbital index r = 
\mj\ = §(§), and 'pseudospin' cr=t(+), !(-)• In this notation, we consider the 'pseudospinor' field 
&(k)=(f^ i , 3 , /fc i _? fl 3 _)' where fl T a is the creation operator for an electron in the 



orbital \mj\=\, § with momentum k and 'pseudospin' a. 
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The representation of the symmetry operations that belongs to the D 4h symmetry of the 
URu 2 Si 2 crystal structure is: T1Z symmetry, inversion symmetry X, four-fold rotational sym- 
metry C 4 , and the two reflection symmetries V x / y . The SO /-state of actinides is invariant under 
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all symmetries except the mirror reflection, which in fact allows the formation of the SO density 
wave into a finite gap in the HO state (see SI). On the basis of these symmetry considerations, it is 
possible to deduce the general form of the non-interacting Hamiltonian as: 

*„-£4t W ( Wfc)M * ) W), (i) 

fe,<7 \ h 2 i(k) h 22 (k) J 

h TT , (k) = e TT , (k)r° + d TT , (k) • r . (2) 

Here, r M (ji e 0,x,y,z) depict the 2D Pauli matrices in the orbital space and r° is the iden- 
tity matrix (a^ matrices will be used later to define the spin space). The T1Z invariance re- 
quires that h 2 2/2i{k)=h\ 1 i 12 {—k). Under T1Z and X, the symmetry of e TT >(k) and d x T ^ z {k) 
must complement to their corresponding identity and Pauli Matrix counterparts, respectively. 
Hence we obtain the Slater-Koster hopping terms as: [e(fe), d x , d y , d z ]u = [ — 2t( cos k x + 
cosky) — /i, — 2tisin/c x , — 2iisin/c y , — 2t 2 (cosk x — cosfc y )], and [e(k), d x , d y , d z \\ 2 = 
[0, 0, 0, — 4:t z cos (k x /2) cos (k y /2) cos (k z /2)] . The obtained values of the tight-binding hopping 
parameters as (i, £ 1? t 2 , ^)=(-45 ,45,50,-25) in meV. The above Hamiltonian can be solved analyti- 
cally which gives rise to four SO-split energy dispersions as 



BT{k) = e(k) + |df 2 (k)P + \d^{k)\\ (3) 

Here a = ± and r — ± become band indices. An important difference of the present Hamiltonian 
with that of bulk topological insulators [0 or quantum spin-Hall systems [3 j is the absence of a mass 
or gap parameter in the former case. The computed non-interacting bands are plotted in Fig.[2fb), 
which exhibit several Dirac points along the high-symmetry lines. Focusing on the Dirac point 
close to Ep 9 we find that it occurs at the crossing between bands E + ~ and E~ + , demonstrating 
that it hosts four-fold Kramer's degeneracy (two orbitals and two spins). Therefore, lifting this 
degeneracy requires the presence of a SO order parameter. However, it is important to note that 
the gap opening at the Dirac point is not a manifestation of the presence of degeneracy at it, but a 
consequence of the SO density wave caused by FS instability. 

SO density wave induced HO. The 'hot-spot' Q h divides the unit cell into a reduced 6 SO 
Brillouin zone' in which we can define the Nambu operator in the usual way (^(fc), &(k + Q%) 9 
W(k + Qh)). In this notation, the modulated SO density wave (SODW) term can be written in 
general as 



Hsodw = : \&(k)T^(k + Q h ) 



(4) 



jJLV 
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where /i, v £ {0, x, y, z}. The symbol :: represents normal ordering. Here g is the contact coupling 
interaction arising from screened interorbital Coulomb term embedded in Hund's coupling param- 
eter, and = r M (g) cr^, r and <r represent Pauli matrices in orbital and spin basis, respectively. 
Absorbing g and T into one term we define the mean-field order parameter 

AP" = g» u (k) ($fi(k)[T t * <g> <r"]V(k + Q h )} . (5) 
= <T(*0 (/^[< T ,0]/ fe+ g,, T/ ^) . (6) 

Here r, r' and a, a' (not in bold font) are the components of the r M and a v matrices, respectively. 
Without any loss of generality we fix the spin orientation along z-directions {y = z). Therefore, we 
drop the index v henceforth. Furthermore we define the gap vector as bP T ,{k) = 5 ,/i A^ r/ (fc)r^ r ,, 
where we split the interaction term g(k) into a constant onsite term and the dimensionless order 
parameter A(fe). With these substitutions, we obtain the final result for the order parameter as 

M " = ( E /jU KAk)<rU] fk+Q h ,r',a\ ■ (7) 

\ttW / 

Eq. [7] admits a plethora of order parameters related to the SO density wave formations which 
break symmetry in different ways. Among them, we rule out those parameters which render gap- 
less states by using the symmetry arguments (see SI): All four order parameters obey 1 symmetry, 
while only M y term is even under T1Z, because it is the product of two odd terms r y and a (we 
drop the superscript V henceforth). This is the only term which commences a finite gap opening 
if the translational or rotational symmetry is spontaneously broken. We have shown in SI that there 
exists a considerably large parameter space of the coupling constant 'g 9 of the coupling constant 
where this order parameter dominates. 

Eq. [7] implies that spin and orbital orderings occur simultaneously along the 4 hot-spot' di- 
rection Q h , as illustrated in Fig. [2^. It propagates along Q x h or Q y h directions with alternating 
signs (particle-hole pairs) to commence a SO density wave. The resulting Hamiltonian breaks the 
four-fold rotational symmetry down to a two-fold one C 2 , and gives rise to a so-called spin-orbit 
'smectic' state which breaks both translational and C 4 symmetry. [45] The present T1Z invariant 
SO order parameter is inherently distinct from any spin or orbital or even interorbital spin-density 
wave order which break T1Z symmetry. This criterion also rules out any similarly between our 
present spin-orbit smectic state with the spin-nematic phase P5ll or spin-liquid state ll36l . Further- 
more, the present order parameter is different from T1Z invariant 'hybridization wave' (between 
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/ and d orbitals of same spin), or charge density wave or others, [21 [71 as SO order involves flip- 
ping of both orbital (between split / orbitals that belong to T 6 symmetry) and spin simultaneously. 
Taking into account the band- structure information that Q h represents the interband nesting, it is 
instructive to focus on only bi 2 (k) component (thus the subscript ' 12' is eliminated hereafter). 
Therefore, the SO density wave does not introduce a spin or orbital moment, but a polarization in 
the total angular momentum Smj=±2 [for the ordering between §(-§) and -|(|)]. 

The b vector belongs to the same irreducible point group representation, E g , of the crystal 
with odd parity, and can be defined by |6(fe) \=2igA x sin k x a, or 2igA y sin k y a for the wavevectors 
Ql = (1± 0.4, 0, 0), or Ql = (0, 1 ± 0.4, 0), respectively. The mean-field Hamiltonian for the HO 
state within an effective two band model reduces to the general form H mf =H +H S odw, where 
the particle-hole coupling term is 

Hsodw = 2i ^2 A ^ sin ( fc M a )/ii,t/*+Q£A+ + h - c - (8) 

k fJ>=x,y 

In the Nambu representation, it is obvious that the HO term merely adds a mass term to the d\ 2 term 
defined above. At the band-crossing points located where \d\ 2 |-|dii|=0, a gap opens by the value 
of |6(fc)| 2 . Figure |2] demonstrates the development of the quasi-particle structure in the HO state. 
The band progression and the associated gap opening is fully consistent with the angle-resolved 
photoemission spectroscopy (ARPES) observations [fT51 [T6l . The scanning tunneling microscopy 
and spectroscopic fSTM/S)ll9l[T4l fingerprints of the gap opening in the density of state (DOS) is 
also described nicely within our calculations, see Fig.[2]c. 



DISCUSSION 



The SO moment is J* rW (q, T) = E fc flrAV[ T rr^aa'}fk+ q y^(T), where T is imaginary 
time. Introducing simplified indices a, (3 = tt'oo', the correlation function of vector can 
be defined as T) = ^ (Tj-J^q, T)Jp(—q, 0)), where Tj- is normal time-ordering. Our 

numerical calculation of the X^piQhiw) within random-phase approximation (RPA) yields an 
inelastic neutron scattering (INS) mode with enhanced intensity at Q h near ujq ~ 4.7 meV below 
T h as shown in Fig. Id. INS data (symbols) at a slightly large momentum agrees well with our 
calculation, however, a polarized INS measurement will be of considerable value to distinguish 
our proposed T1Z invariant mode from any spin-flip and elastic background. [46J A T-dependent 
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study of the INS mode also reveals that this mode becomes strongly enhanced at Qh rather than at 
the commensurate one below TL lfT2l 

One way to characterize the nature of a phase transition is to determine the temperature evolu- 
tion of the gap value. Our computed self-consistent values of the mean-field gap A(T) agree well 
with the extracted gap values from the STM spectra[9j [see inset to Fig. 2c]. In general, the entropy 
loss at a mean-field transition is given by ll20ll AS ~ k B j^, where A is the HO gap and is the 

s F 

Fermi energy of the gapped state. At HO the Fermi energy £ F « yj \ [Ei(k F ) — E 2 (k F )] 2 + A 2 , 
where the two linearly dispersive bands near the Fermi level yields Ei(k F ) w hv iF k iF . Using 
the measured Sommerfield coefficient 7=180 mJmol _1 K -2 , compared to its linear expansion of 
7o=50 mJmol _1 K~ 2 , we obtain the mass renormalization factor Z~ x = 7/70 = 3.6. This gives 
v F = Zv l F . For the two bands that participate in the HO gap opening, we get vi F =0.lla/h in 
eV at kiF=Q.5ir/a and v 2F =0.58a/h in eV at k 2F =03n/a from Fig. la. Using the experimental 
value of A=5meV.ll9in~4l we obtain AS ~ 0.28/c s In 2, which is close the experimental value of 
0.3fc B ln2.ni 

We now evaluate the topological invariant index of interacting Hamiltonian in Eq.[8]to demon- 
strate that HO gap opening in URu 2 Si 2 also induces topological phase transition. To characterize 
the topological phenomena, we recall the Fu-Kane classification scheme [0 which implies that if a 
time-reversal invariant system possess an odd value of Z 2 invariant index, the system is guaranteed 
to be topologically non-trivial. Z 2 index is evaluated by the time-reversal invariant index v { = ±1, 
if defined, for all filled bands as Z 2 = v\v 2 ...v n , where n is the total number of orbitals in the 
Fermi sea. A more efficient method of determining the topological phase is called the adiabatic 
transformation scheme used earlier in realizing a large class of topological systems, especially 
when Z 2 calculation is difficult. [47 j In this method, the non-trivial topological phase of a system 
can be realized by comparing its band-progression with respect to an equivalent trivial topological 
system. URu 2 Si 2 is topologically trivial above the HO state, i.e. Z\ = +1. The gap opening 
makes the top of the valence band (odd parity) to drop below E F as shown in Fig. [2]b. Thereby, 
an odd parity gained in the occupied level endows the system to a non-trivial topological metal. 
To see that we evaluate the topological index for the HO term as v ho = J dfcfi(fe), where the 
corresponding Berry curvature can be written in terms of 6- vector as Q = b • x for each 
spin with b = b/\b\. Due to the odd parity symmetry of 6, it is easy to show that v ho — — 1 which 
makes the total Z 2 value of the HO phase to be Z 2 = Z 2 ° x ^ = -1, and hence we show that 
hidden order gapping is a topologically non-trivial phase. The consequence of a topological bulk 



9 



gap is the presence of surface states. 0147]] In our present model, we expect two surface states of 
opposite spin connecting different orbitals inside the HO gap. As the system is a weak-topological 
system, the surface states are unlikely to be topologically protected. The spin-orbit locking of 
these states can be probed by ARPES using circular polarized incident photon which will be a 
definite test of this postulate. 

The HO gap is protected from any T1Z invariant perturbation such as pressure (with suffi- 
cient pressure the HO transforms into the LMAF phase), while T1Z breaking perturbation such 
as magnetic field will destroy the order. Remarkably, these are the hallmark features of the HO 
states, HH [23l [35) which find a natural explanation within our SO density wave order scenario. 
In what follows, the magnetic field will destroy the HO state even at T = K, that means at a 
quantum critical point (QCP) as the HO is a spontaneously broken symmetry phase. lfT9ll However, 
due to the finite gap opening at the HO state, it requires finite field to destroy the order. The ther- 
modynamical critical field can be obtained from[Q~5]] A = (xg(^ re5 ))S^, where B c is the critical 
field and XQ^res) ~ 2Aa 2 tanh(A/2k B T h )/cc^ es at the resonance mode that develops in the HO 
state, a = g^B\(dmj)\ = 2g/ie and the bare 4 g'-factor g = 0.8. Substituting uo res = 4.7meV, we 
get the location of the QCP at B w 38 T, which is close the experimental value of B = 34 T. [fT8ll 

Broken symmetry FS reconstruction leads to enhanced Nernst signal. [5 lj For the case of bro- 
ken symmetry spin-orbit order, we expect to generate spin-resolved Nernst effect which can be 
measured in future experiments to verify our proposal. fl52) 

In summary, we proposed a novel SO density wave order parameter for the HO state in 
URu 2 Si 2 . Such order parameter is T1Z symmetry invariant. We find no fundamental reason 
why such order parameter cannot develop in other systems in which both electronic correlation 
and SO of any kind are strong. Some of the possible materials include heavy fermion systems, 
Iridates,L53j SrTi0 3 surface states, fl54) SrTi0 3 /LiA10 3 interface, fl55) Half-Heusler topological 
insulator ll47lh and other d- and /-electron systems with strong SO. In particular, a Rashba-type 
SO appears due to relativistic effect in two-dimensional electron system yielding helical FSs. In 
such systems, the FS instability may render similar SO density wave, and the resulting quasiparti- 
cle gap opening is observed on the surface state of BiAg 2 alloys even when the spin-degeneracy 
remains intact. Il56ll Furthermore, recent experimental findings of quasiparticle gapping in the 
surface state of topological insulator due to quantum phase transition even in the absence of time- 
reversal symmetry breaking can also be interpreted as the development of some sort of spin orbit 
order. 11571 



10 



[1] Bernevig, B. A., Hughes, T. L., Zhang, S. C. Quantum spin Hall effect and topological phase transition 

in HgTe quantum wells. Science 314, 1757 (2006). 
[2] Fu, L., Kane, C. L., Mele, E. J. Topological insulators in three dimensions. Phys. Rev. Lett. 98, 106803 

(2007) . 

[3] Zhang H. et al Topological insulators in Bi2Se3, Bi2Te3 and Sb2Te3 with a single Dirac cone on the 

surface. Nat. Phys. 5, 438 (2009). 
[4] Hsieh, D. et al. A tunable topological insulator in the spin helical Dirac transport regime. Nature 460, 

1101 (2009). 

[5] Raghu S, Qi XL, Honerkamp C, Zhang SC Topological Mott insulators. Phys. Rev. Lett. 100, 156401 

(2008) . 

[6] Dzero, M., Sun, K., Galitski, V., Coleman, P. Topological Kondo insulators. Phys. Rev. Lett. 104, 
106408 (2010). 

[7] Mong, R. S. K., Essin, A. M., Moore, J. E. Antiferromagnetic topological insulators. Phys. Rev. B 81, 
245209 (2010). 

[8] Mydosh, J. A., Oppeneer, R M. Colloquium: Hidden order, superconductivity and magnetism - the 

unsolved case of URu 2 Si 2 . Rev. Mod. Phys. 83, 1301 (2011). 
[9] Aynajiana, R et al. Visualizing the formation of the Kondo lattice and the hidden order in URu2Si2. 

Proc. Natl. Acad. Sci. USA 107,10383 (2010). 
[10] Palstra, T. T. M. et al. Superconducting and magnetic transitions in the heavy-fermion system 

URu 2 Si 2 . Phys. Rev. Lett. 55, 2727 (1985). 
[11] Palstra, T. T. M. et al, Anisotropic electrical resistivity of the magnetic heavy-fermion superconductor 

URu 2 Si 2 . Phys. Rev. B 33, 6527 (1986). 
[12] Ramirez, A. P. et al. Nonlinear susceptibility as a probe of tensor spin order in URu 2 Si 2 . Phys. Rev. 

Lett. 68, 2680 (1992). 

[13] Visser, A. D. et al. Thermal expansion and specific heat of monocrystalline URu 2 Si 2 . Phys. Rev. B 34, 
8168 (1986). 

[14] Schmidt, A. R. et al. Imaging the Fano lattice to a hidden order transition in URu 2 Si 2 . Nature 465, 
570 (2010). 

[15] Santander-Syro, A. F. et al. Fermi-surface instability at the 'hidden-order' transition of URu 2 Si 2 . Nat. 



11 



Phys. 5, 637 (2009). 

[16] Yoshida, R. et al. Signature of hidden order and evidence for periodicity modification in URii2Si2. 

Phys. Rev. B 82, 205108 (2010). 
[17] Mason, T. E., Buyers, W. J. L. Spin excitations and the electronic-specific heat of URu2Si2. Phys. Rev. 

£43, 11471 (1991). 

[18] Broholm, C. et al. Magnetic excitations and ordering in the heavy-electron superconductor URu2Si2. 

Phys. Rev. B 43, 12809(1991). 
[19] Mentink, S. A. M. et al. Gap formation and magnetic ordering in URu2Si2 probed by high- field mag- 

netoresistance. Phys. Rev. B 53, 6014(R) (1996). 
[20] Chandra, R et al. Hidden orbital order in the heavy fermion metal URu2Si2. Nature 417, 831-834 

(2002). 

[21] Amitsuka, H. et al. Effect of pressure on tiny antiferromagnetic moment in the heavy-electron com- 
pound URu 2 Si 2 . Phys. Rev. Lett. 83, 5114 (1999). 

[22] Fisher, R. A. et al. Specific heat of URu2Si2i Effect of pressure and magnetic field on the magnetic 
and superconducting transitions Physica B 163, 419 (1990). 

[23] Okazaki, R. et al. Rotational symmetry breaking in the hidden-order phase of URu2Si2. Science 331, 
439 (2011). 

[24] Barzykin, V., Gofkov, L. R Possibility of observation of nontrivial magnetic order by elastic neutron 

scattering in magnetic field. Phys. Rev. Lett. 70, 2479 (1993). 
[25] Kasuya, T. et al. Hidden ordering and heavy mass in URu2Si2 and its alloys. /. Phys. Soc. Jpn. 66, 

3348 (1997). 

[26] Ikeda, H., Ohashi, Y. Theory of unconventional spin density wave: A possible mechanism of the 
micromagnetism in U-based heavy fermion compounds. Phys. Rev. Lett. 81, 3723 (1998). 

[27] Mineev, V.R, Zhitomirsky, M. E. Interplay between spin-density wave and induced local moments in 
URu 2 Si 2 . Phys. Rev. B 72, 014432 (2005). 

[28] Okuno, Y, Miyake, K. Induced-moment weak antiferromagnetism and orbital order on the itinerant- 
localized duality model with nested FS : A possible origin of exotic magnetism in URu2Si2. /. Phys. 
Soc. Jpn. 67, 2469 (1998). 

[29] Varma, C. M., Zhu, L. Helicity order: Hidden order parameter in URu 2 Si 2 . Phys. Rev. Lett. 96, 036405 
(2006). 

[30] Santini, P. Behavior of URu 2 Si 2 in an applied magnetic field. Phys. Rev. B 57, 5191-5199 (1998). 



12 



[31] Hanzawa, K. Hidden octupole order in URu 2 Si 2 . /. Phys.: Condens. Matter 19, 072202 (2007). 
[32] Haule, K., Kotliar, G. Arrested Kondo effect and hidden order in URu 2 Si 2 . Nat. Phys. 5, 796-799 
(2009). 

[33] Elgazzar, S. et al. Hidden order in URu 2 Si 2 originates from FS gapping induced by dynamic symmetry 

breaking. Nat. Mat. 8, 337-341 (2009). 
[34] Dubi, Y., Balatsky, A. V. Hybridization wave as the hidden order in URu 2 Si 2 . Phys. Rev. Lett. 106, 

086401 (2011). 

[35] Fujimoto, S. Spin nematic state as a candidate of the hidden order phase of URu 2 Si 2 . Phys. Rev. Lett. 
106, 196407 (2011). 

[36] Pepin, C. et al. Modulated spin liquid: A new paradigm for URu 2 Si 2 . Phys. Rev. Lett. 106, 106601 
(2011). 

[37] Oppeneer, R M. et al. Spin and orbital hybridization at specifically nested FSs in URu 2 Si 2 . Phys. Rev. 
£84, 241102(R) (2011). 

[38] Riseborough, R S. et al. Phase transition arising from the underscreened Anderson lattice model: A 
candidate concept for explaining hidden order in URu 2 Si 2 . Phys. Rev. B 85, 1651 16 (2012). 

[39] Wiebe, C. R. et al. Gapped itinerant spin excitations account for missing entropy in the hidden-order 
state of URu 2 Si 2 . Nat. Phys. 3, 96 (2007). 

[40] Blaha, P. et al. An augmented plane wave plus local orbitals program for calculating crystal properties. 
(Techn. Univ. Wien) (2001). 

[41] Perdew, J. P. et al. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865 (1996). 

[42] Denlinger, J. D. et al. Comparative Study of the Electronic Structure of XRu 2 Si 2 : Probing the Ander- 
son Lattice /. Elect. Spect. Relat. Phenom. 117, 347 (2001). 

[43] Hotta, T., Ueda, K. Construction of a microscopic model for / electron systems on the basis of a j — j 
coupling scheme Phys. Rev. B 67, 104518 (2003). 

[44] Winkler, R. Spin-orbit coupling effects in two dimensional electron and hole systems. (Springer Tracts 
in Modern Physics, Vol. 191, Springer (2003). 

[45] Zaanen, J., Gunnarsson, O. Charged magnetic domain lines and the magnetism of high-T c oxides. 
Phys. Rev. B 40, 7391-394 (1989). 

[46] Fauque, B. et al. Magnetic Order in the Pseudogap Phase of High-TC Superconductors. Phys. Rev. 
Lett. 96, 197001 (2006). 

[47] Lin, H. et al. Half-Heusler ternary compounds as new multifunctional experimental platforms for 



13 



topological quantum phenomena. Nature Materials 9, 546 (2010). 
[48] Kim, K. H. et al. Magnetic-Field-Induced Quantum Critical Point and Competing Order Parameters 

in URu 2 Si 2 . Phys. Rev. Lett. 91, 256401 (2003). 
[49] Sachdev. S. Quantum phase transitions, Cambridge University Press, Cambridge U.K. (1999). 
[50] Pathria, R. K. Statistical Mechanics, Butterworth-Heinemann publishers, 2nd Edition, (1996). 
[51] Oganesyan, V. & Ussishkin, I., Nernst effect, quasiparticles, and d— density waves in cuprates. Phys. 

Rev. B 70, 054503 (2004). 

[52] Cheng, S.-G., Xing, Y., & Sun, Q. & Xie, X. C. Spin Nernst effect and Nernst effect in two- 
dimensional electron systems. Phys. Rev. B 78, 045302 (2008). 

[53] Witczak-Krempa, W., Kim, Y.B. Topological and magnetic phases of interacting electrons in the py- 
rochlore iridates. Phys. Rev. B 85, 045124 (2012). 

[54] Santander-Syro, A. F. et al. Two-dimensional electron gas with universal subbands at the surface of 
SrTi0 3 . Nature 469, 189-193 (2011). 

[55] Caviglia, A. D. et al. Tunable Rashba spin-rbit interaction at Oxide interfaces. Phys. Rev. Lett. 104, 
126803 (2010). 

[56] Bentmann, H. et al. Direct observation of interband spin-orbit coupling in a two-dimensional electron 

system. Phys. Rev. Lett. 108, 196801 (2012). 
[57] Sato, E. et al. Unexpected mass acquisition of Dirac fermions at the quantum phase transition of a 

topological insulator. Nat. Phys. 7, 840 (2011). 
[58] Haraldsen, J. T. et al. Hidden-order pseudogap in URu 2 Si 2 . Phys. Rev. B 84, 214410 (201 1). 
[59] Bourdarot, F. et al Precise Study of the Resonance at Qq = (1, 0, 0) in URu2Si2. /. Phys. Soc. Jap. 

79, 064719(2010). 
Acknowledgments 

The author thanks A. V. Balatsky, M. J. Graf, A. Bansil, R. S. Markiewicz, T. Durakiewicz, J.-X. 
Zhu, P. M. Oppeneer, J. Mydosh, P. Wolfle, and J. Haraldsen for useful discussions. Work at the 
Los Alamos National Laboratory was supported by the U.S. DOE under contract no. DE-AC52- 
06NA25396 through the Office of Basic Energy Sciences and the UC Lab Research Program, and 
benefited from the allocation of supercomputer time at NERSC. 



14 




FIG. 1. Ab initio band structure and Fermi surface of URu2Si2. a, Computed non-interacting energy 
dispersions of URu 2 Si 2 , using Wien2K software, fiOlHTll are presented along T(0,0,0), X(7r,0,0), M(7r,7r,0), 
and Z(0,0,7r) directions. The band structure is consistent with the previous full potential local orbitals 
(FPLO) and full potential linearized augmented plane wave (FPLAPW) calculations in the paramagnetic 
state. 1 33 1 The low-energy dispersions along T-X is expanded in b and contrasted with the same but without 
the SO coupling in c. The FS in the k z = plane is shown in d. The red arrow dictates the FS 'hot-spot' 
that emerges after including SO coupling. 
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FIG. 2. Spin-orbit density wave and the hidden-order gap opening, a, A typical form of the staggered 
SO order is schematically described for an illustrative case of commensurate wavevector. The solid and 
dashed circles encode two opposite orbitals, r = ±, where the associated arrows depict their 'pseudospins' 
a. Both r and cr, representing orbital and spin respectively, individually break T1Z symmetry, while their 
product remains T1Z invariant, b, Model dispersions of the | ± ^), and | ± |) subbands plotted along the 
axial direction. Black and red lines give dispersion before and after including the HO gap, respectively. 
An artificially large value of A=50 meV is chosen here to clearly explicate the momentum dependence of 
the modulated SO gap opening, c, Modifications of DOS upon entering into the HO phase are compared 
with measured DOS in the STM experiment (green line). 03) Note that the experimental data is subtracted 
from the background spectrum at T>T^, which helps highlight the appearance of multiple structures in 
the DOS spectrum at the HO state. Here the gap magnitude A(0)=5 meV, obtained at a coupling strength 
of g = 27 meV, see SI. Inset: The self-consistent value of A(T) exhibits the mean-field behavior of the 
HO gap, in consistent with experiments. ||9l We obtain = 22 K which is larger than the experimental 
value of Th = 17.5 K. However, recently it has been pointed out that there exists a 'pseudogap' above the 
HO state,[ 17 1 which presumably reduces the mean-field temperature scale, d, RPA result of SO correlation 
function at g =28.4 meV shows a resonance peak at ujq = 4.7 meV at Qh, in good agreement with 
experimental data. ifTUl fTTTl 



16 



SUPPLEMENTARY INFORMATION (SI) FOR "SPIN-ORBIT DENSITY WAVE INDUCED 
HIDDEN TOPOLOGICAL ORDER IN URU 2 SI 2 " 

In the main text, we have provided an effective two band model which is relevant for the study 
of hidden-order (HO) gap structure coming from Fermi surface (FS) instability. Here, we expand 
on how such an effective Hamiltonian is deduced. The spin-orbit splitting state is, in general, stud- 
ied within either L.S coupling or j — j coupling approaches. In the former case, the total L and 
S are formed due to strong Hund's coupling prior to the formation of spin-orbit coupled eigen- 
states. Such process mainly occurs in insulating compounds with localized / states. [lj However, 
in actinides the spin-orbit (SO) coupling is stronger than the Hund's coupling. [2] Therefore, total 
angular momentum J = L + S is the good quantum number for this state. In such systems, a 
SO density wave order in the two-particle channel can arise at some critical value of the coupling 
constant g, by taking advantage of any instability, such as FS 'nesting' (shown in the main text), 
even when the time-reversal symmetry remains invariant. The coupling g can be related to some 
form of 'screened' interorbital Coulomb term. Physically, SO density wave is different from a 
spin density wave because here spin-flip occurs between two different orbitals without breaking 
time-reversal symmetry. Therefore, we desire to study a Hamiltonian: 

H = Hq + Hsodw, (9) 

where H is the non-interacting part and H S odw is the SO density wave term. In the 'pseudospin' 
basis introduced in the main text &(k) = i (7 , f ] k 3 ^ f ] k i ^ f ] k z -) (a = -o), the repre- 
sentation of the symmetry operations for URu 2 Si 2 system which belong to the D 4h symmetry is: 
time-reversal symmetry T1Z = K • icr y ® r°, inversion symmetry 1 = r z ® r°, four-fold ro- 
tation symmetry around the z axis C 4 = exp [i(n/4:)a z ® r°] and the two reflection symmetries 
V x / y = icr x / y which map x — )► — x (where x is in T-X direction) and y — >• —y (where y is or- 
thogonal to T-X direction), respectively. Here, /C is complex conjugation operator, and a x,y,z and 
T x, y ,z ^p^t the two-dimensional Pauli matrices in the 'pseudospin' and orbital space, respectively 
where r° is the unitary matrix. 

Each symmetry operation transforms the time-reversal invariant /-electron field as: 

X^{k X) ky) -+ -T 4f(-k X , ~ky), (10) 

Cjf(k x ,k y ) ^ir y 4>(-k y ,k x ), (11) 
Vx/ y *(k x , k y ) -> TT z V(Tk x , ±k y ). (12) 
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These symmetry operations imply that the spin-orbit coupled actinide /-state is an odd-parity 
wavefunction. Using these symmetry properties, the non-interacting Hamiltonian is deduced in 
the main text, using standard procedure, see for example Refs.|3]and|4[ 



HIDDEN ORDER PARAMETER 

Considering the non-interacting FS nesting at Q h = (1 ± 0.4, 0, 0) for the hidden order state 
as demonstrated in the main text, we expand the above-mentioned 'pseudospinor' in the Nambu 
notation as &(fe) = (^(fc), &(fe + Qh))- I n this notation, the modulated SOC interaction term 



can be written in general as 

r - , - i 2 

(13) 



Hsodw = YY v - \&(k)T^(k + Q h j\ 



[IV 



where /x, v £ {0, y, z}. The symbol :: represents normal ordering. Here g is the coupling 
constant discussed latter and = r M ® <r^ , r and <r represent Pauli matrices in orbital and spin 
basis, respectively. Absorbing g and T into one term we define the mean-field order parameter 



M w = ^( fc ) /^r(fe)t[ r /x ^(fc + . (14) 

= <rik) (fU a i<> ® w,^) . d5) 



Here r, r' and a, a 7 (not in bold font) are the components of the r M and cr v matrices, respectively. 
Without any loss of generality we fix the spin orientation along z-directions {y — z). Therefore, 
we drop the index v henceforth. Furthermore we define the gap vector as 

6^(fc) = A^(fc)< T „ (16) 

where momentum dependence of the gap function transforms according to the same irreducible 
representation of the point-group symmetry (g is absorbed in the gap function defined below). 
With these substitutions, we obtain the final result for the order parameter as 

^ M = ( E /jU KAk)* z a& ] f k+Qh y,a) • (17) 

\ttW / 

For the unidirectional modulation vector Q h , all components of M M break C 4 rotational symme- 
try. All the symmetry properties of the order parameters are given in Table I. M°, M x and M z 
break time-reversal and thus are ruled out as the hidden-order state is arguably does not exhibit 
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any time-reversal symmetry breaking. [[5]43. Although some evidences for the time-reversal sym- 
metry breaking are also present, but it is not well established if the measurements are done in 
single crystal where the time-reversal symmetry breaking LMAF state and time-reversal symme- 
try invariant HO state are not mixed. In microscopic sense, LMAF and HO state are separated by 
phase transition, and thus both phases cannot inherit same broken symmetry. 





M° 


M x 


M y 


M z 


m 






+ 
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+ 


+ 
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'Px/y 


+ 






+ 



TABLE I. Symmetry properties of various density wave parameters. '+' ('-') represent 'even' ('odd') parity 
of the order parameter under the corresponding symmetry operations. 

On the other hand, M 2 is even under time-reversal but odd under parity, and thus adds a mass 
term to the Hamiltonian which opens a gap at the nested portion of the Fermi surface. This is the 
term that represents modulated spin-orbit coupling. A trivial check can be performed by explicitly 
writing down the the spin-orbit coupling term ^i(-l) ix L iiZ Si iZ oc 4+ ><r /*,+ |0 -' - fi-^k-^' = 
2i f{iji,l*' + h ' c -> where we substituted f ii±i<T = (/^ ± i/^). 

With this definition, the total Hamiltonian in the reduced Brillouin zone can be written as 

H= £ (h TT ,{k)fl T J k , Tl ,-* 

+ h TT ,{k + Qh)fl +Qh , T ,Jk+Q h ,T',a) 

I 

- J2 { M ^)fl,rA< Tl ®a: a ,}h + Q h ,r'y 

+ M»(k + Qh)fk+Q h ,T,a \Ttt' ® <Afw) 

+ ^M^(k)M^(k + Q h )/g^. (18) 

Where h is the non-interacting part of the Hamiltonian defined in Eq. 2 of the main text. The 
prime over a summation implies that the summation is performed in the reduced Brillouin zone. 
Diagonalizing the above Hamiltonian, we deduce the quasiparticle states are: E]? v = + vE]?^, 
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FIG. 3. a, Black and blue symbols are ab-initio FSs, as shown in Fig. 1 of main text, but plotted here only 
two relevant bands. Two large square FSs for band 1 are shifted along the x— direction by nesting vectors 
Qn = Q2 and Q12, while the two small square FSs for band 2 are shifted along — x direction. Good 
nestings for all three vectors are observed, but the intraband one are protected by strong SOC, where Q12 = 
Qh is the hidden order nesting that opens a gap. b, Phase diagram as a function of intra-orbital {g^j = 
9012 = 5021) an d inter-orbital coupling constant (gou = #011 = #022) and nearest neighbor hopping (t) 
with respect to the Fermi energy (Ep) for all thee nesting instabilities. We find that for smaller inter-orbital 
coupling than the intra-orbital one, the hidden order parameter arises if Epjt ration is large. With large t 
(proportional to pressure), the LMAF phase at the commensurate nesting wins. This result is consistent, at 
least, in qualitative level, with the phase diagram of URu^Si^.fHl 



and E\ 



V(a-) 2 + |A.r'| 2 , where = ± £^ Q )/2, and 



E r *(k) = e(k)+r /^K 2 (k)|2 + (7 



Ei<i( fe )i 2 



(19) 



from Eq. 3 of the main text, u, v' = ±. The corresponding coherence factors are 




(20) 
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SELF-CONSISTENT GAP EQUATION 



We can easily derive the self-consistent gap equation from Eq.[T7]using Bogolyubov treatment. 
We substitute the fermion operator / in terms of a Bogolyubov operators as 

/TcT TcT 
k,ra = U k 7fc 5T(7 - V k 7fc+Q hjT0 - 

fk+Q h ,ra = ^klk.ra + U™'Yk+Q hi T<T' (21) 

From Eqs.[l5j[T6} andfTTJ we get the gap function as 

Kh = C> E (ft +Qh ,rJ^) , (22) 

k,a 



k,a 

" l k u k 



ul'vT G 



k.a 



^k+Q h ,T^k+Q h ya 
k,a 



k.a 



J k 



(23) 



Here we have substituted (7^ Ta r ) k ^ a ) = (fe), where n is the Fermi function. We find that 
#12 = 27 meV gives A =5 meV in consistent with experiments. For this temperature independent 
value of ^12, we obtain T h = 22K, which is higher than the experimental value of T h = 17.5 meV. 
Possible reason for overestimating the value of T h are the neglect of quantum fluctuation, disorder 
which can reduce its value. In this context, it can be noted that recently, a 'pseudogap' phase upto 
20 K is marked from experimental features, Q3 which arguably suggests that there are indeed 
fluctuations present above T h . 

SPIN-ORBIT CORRELATION FUNCTION 

The non-interacting single-particle Green's function for the Hamiltonian given in Eq. 1 above is 
defined as G rrW (fc, iuj n ) = -1/(3 T,n( T rfkra(T)fl r ; a ;(0))e iu;nT . Here T is the imaginary time, 
f kra (T) is the imaginary time evolution of the fermionic operator f kra , n is the fermionic Matsub- 
ara frequency, and f3 = l/k B T, where k B is Boltzmann constant. Tj- gives normal time-ordering. 
The anomalous part ofthe Green's function is F TT , aa ,{k,iu n ) = -J ^ n \Trfkra(T)fl k+Q)r/(7/ (0) 
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By Fourier transforming in the Matsubara frequency space, we obtain the explicit forms of two 
Green's functions as 



G TT ' aa '(k,iuj n ) = £i^ rW ,(fc) 



^ (24, 



iu n - El* iu n - El* 



F Tr , aa ,(k,iu; n ) = ^i?^ w (fc)^X V 

1 1 



iu n - E™ iu n - El* 



(25) 



The orbital overlap matrix-element R v TT , a(J ,{k) = ^ cr (fc)^j cr ,(fc), where ip is the eigenvector of 
the noninteracting Hamiltonian, projects the Green's function from band basis to the orbital one. 
However, to simplify our calculation, we assume that each block state E T * corresponds to each 
orbital which makes R = 1. By taking analytical continuation of the Matsubara frequency to the 
real frequency in the above Green's function, it is easy to show that the gap function A can be 
obtained by averaging the anomalous Green's function. 
The general form of the polarization vector is 

Jrr'.AQ, T) = £ ftrAT)T» r ,*l a ,f k+q , Tl , a ,{T). (26) 
k 

To simplify the notations, we define composite indices a, (3 = tt'oo' , in which the correlation 
function of vector can now be defined as 

&q,r) = ± (T T J z a (q,T)J z e(-q,0)) . (27) 

Substituting J z a and then applying standard Wick's decomposition to the electron brackets yields, 

Xap(q, Wm) = -jj £ [G a (k, iu n )G p {k + q, iu n + ip m ) 

k,n 

-F a (k, ioj n )F p (k + q, iuj n + ip m ))\. (28) 



= 4£ 5 ^(fc,q) 



N 

k 



x [{u k u k+q + v k v k+q ) (x ++ + X ) 
+ {u k v k+q - v k u kq ) 2 (x + ~ + X~ + ) 



(29) 



(30) 



Here S a p(k, q) = R a (k)Rp(k + q) is the matrix-element term which projects the susceptibility 
from the band representation to the orbital one. x vv ' are the fermionic oscillator terms in the band 



space v, v 



' = ±: 




nK(k) -nK'(k + q) 



(31) 



X 



The momentum, energy, a indices on both sides of the above equation are implicit. 



Hidden order instability 



The divergence in the real-part of the Xa/3(q,u = 0) indicates an instability due to the FS 
nestings between two orbitals a and /3 (in the explicit notation xn means Xim> etc )- The two 
intraband nestings at Q n and Q22 can give rise to spin-density wave or antiferromagnetism, if the 
time-reversal symmetry is broken. We believe that Qn is responsible for the LMAF phase. Our 
present interest is the interband one at Q12, shown in supplementary Fig. [3jA) which happens at 
(1 ± 0.4, 0, 0). Since in this case, both orbital and spin flip together, the time-reversal symmetry 
remains intact. 

The corresponding critical value of the coupling constants g a p at which a gap opening or an 
order parameter develops can be evaluated within random-phase approximation (RPA). The stoner 
criterion for an instable state implies that 1 — XapiQapifydap > 0, or g a/3 < 1 /x'apiQ a/3,0) • 
Looking at the FS areas for each band, shown in supplementary Fig. [3]a, we immediately see 
that l/xii(Qn, 0) < 1/X2 2 (Q22, 0) < l/xi 2 (Qi2j 0), leading to a phase diagram shown in sup- 
plementary Fig. ^B). The present calculation does not incorporate the possible coexistence state 
between different phases. The phase diagram implies that there is a considerably large parameter 
space, where the Q i2 = Qh nesting dominates. We have not considered all possible phases except 
M y for Q u nestings discussed in Eq.[T7] above, because constrained by the symmetry arguments 
given in Table I, others render gapless state. 



scattering mode with enhanced intensity at Q h whose energy scale is approximately given by 



x"(Q/i, u) - CS(uj -El- E 2 k+Q ) w C6(uj - \A(k)\ - \A(k + Q h )\). C is the prefactor which 



has to be evaluated rigorously, but it does not contribute to the peak position in bare x" Here 
we have substituted the condition that non-interacting bands are nested on the FS at Q h such as 



Neutron mode at the Hidden order state 




accompanies an inelastic neutron 
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FIG. 4. The real and imaginary part of the bare susceptibility calculated from Eqs. [30pT] in the hidden 
order state. The red line gives the INS mode obtained using RPA calculation at gu = 28.4 meV whose 
inverse is shown by horizontal dashed line. The RPA value has a huge peak (resonance) at uoq =4.7 meV 
when x'iQh, wq) = l/gu, and we have rescaled its intensity arbitrarily to fit it into the same figure with 
the bare values for comparison. The obtained resonance is in good agreement with experiment shown by 
symbols. ifTTTl 

= |A(fc)| ~5 meV. The RPA correction shifts the resonance energy which depends on the 
value of coupling constant. At gr 12 = 28.4 meV, the l/gu line cuts twice to x'iQh, w) 9 however, 
the resonance is stronger at the energy where the broadening function coming from x" is weaker. 
Therefore, the strong intensity or a resonance occurs at ujq ~ 4.7 meV, see Fig. |4| INS measure- 
ments have observed this resonance at the incommensurate nesting vector Q = (1.4, 0, 0). lfTlT - [T3ll 
However, one have to be careful to directly compare our result with this data. Because, in the 
present case, we expect a mode which does not break time-reversal symmetry. Therefore, in order 
to observe this mode, one requires to perform a polarized INS measurement. 

On the same reasoning, INS should see more resonances at Q n and Q22 even in the non- 
polarized condition. However, the energy scale and intensity of those modes will depend on the 
location of the phase diagram of URu 2 Si 2 where the experiment is performed. As we argued 
earlier, Qn is most likely responsible for the LMAF phase, therefore, we can expect a mode at 
twice of the corresponding gap opening. The peak in the INS spectra at the commensurate vector 
is also observed at a much lower in energy, however, the peak is much broader than the resonance 
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FIG. 5. Proposed phase diagram of URii2Si2 along the magnetic field (B) direction. The data along 
pressure (P) — T plane is taken from Ref . [U The symbols along B — T axes are extracted from Ref . [TU A 
QCP at T = along the field directions is expected from our theory, and also observed in experiment. HH 
Deducing the phase diagram for SC and other possible phases that may arise above the QCP is beyond the 
scope of the present calculation. 



peak observed at the incommensurate vector. lfTTl[T3ll 



QUANTUM CRITICAL POINT 



As mentioned in the main text, any TTZ breaking perturbation such as magnetic field will de- 
stroy the TTZ invariant HO phase. As the HO state incipiently is a spontaneously broken symmetry, 
it will exhibit a second order phase transition along the field axis. In what follows we expect to 
obtain a quantum critical point (QCP) at T = 0, as extensively proposed to be associated with 
any second order phase transition, lfT9l along the magnetic field axis. We find the critical field to 
be B ^38 T, whereas the experimental data lfT8l indicates that the QCP resides around B ~ 34 T. 
The present model can not deduce the phase diagram for the superconducting (SC) state, possibly 
intervening the HO state, or any other phases that may arise above the QCP. [[18) More experimen- 
tal data and theoretical modelling are required to understand the details of this regime of the phase 
diagram. 
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